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Liquid metal embrittlement (LME) is a common feature of systems in which a low melting point 
liquid metal is in contact with another, higher melting point, poly crystalline metal. While different 
systems exhibit different LME fracture characteristics, the penetration of nanometer-thick liquid 
metal films along the grain boundary is one of the hallmarks of the process. We employ EAM 
potentials optimized for Al-Ga binary alloys in a series of MD simulations of an Al bicrystal (with 
a S5 36.9°/[010] symmetric tilt boundary) in contact with liquid Ga with and without an applied 
stress. Our simulations clarify the mechanism of LME and how it is affected by applied stresses. 
The interplay of stress and penetrating Ga atoms leads to the nucleation of a train of dislocations 
on the grain boundary below the liquid groove root which climbs down the grain boundary at a 
nearly constant rate. The dislocation climb mechanism and the Ga penetration are coupled. While 
the dislocations do relax part of the applied stress, the residual stresses keep the grain boundary 
open, thereby allowing more, fast Ga transport to the penetration front (i.e., Ga layer thickening 
process). The coupled Ga transport and "dislocation climb" is the key to the anomalously fast, 
time-independent penetration of Ga along grain boundaries in Al. The simulations explain a wide 
range of experimental observations of LME in Al-Ga the literature. 

PACS numbers: 62.20.Mk, 68.08.De, 81.40.Np 



G 

O 

o 



> 

o 

m 
m 
O 

o 

-i— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

There are many examples in which very deep grooves 
form at the intersections of grain boundaries and the sur- 
face in systems in which a liquid metal is in contact with 
a polycrystalline solid. In some systems, such as Al-Ga, 
Zn-Ga, Cu-Bi and Ni-Bi, the liquid film quickly pene- 
trates deep into the solid along the grain boundary and 
leads to brittle intergranular fracture under the influence 
of even modest stresses. This is a form of liquid metal em- 
brittlement (LME). Although this phenomenon is quite 
common in material processing, LME is not well under- 
stood. LME is particularly important in nuclear reactor 
scenarios in which liquid metals are used as coolants and 
as spallation targets. 

The Al-Ga couple is a particularly well-known LME 
system. Several studies have shown that the maximum 
load a polycrystalline Al sample in contact with liquid 
Ga can sustain decreases as the quantity of Ga on the 
grain boundaries increases (characterized by exposure 
time of these Al samples to Ga) eventually leading to 
intergranular brittle fracture. 1-3 Transmission electron 
microscopy 4-6 (TEM), scanning electron microscopy 7,8 
(SEM), and synchrotron radiation micro-radiography 
studies 9-12 all show that liquid Ga penetrates into grain 
boundaries in Al at a remarkable rate, leading to a dis- 
tinct channel morphology. The penetration of liquid Ga 
along the grain boundaries produces wetting layers with 
thickness ranging from several monolayers 4-6 to several 
hundred nanometers, 10-12 even in the absence of an ap- 
plied load. Interestingly, the rate of propagation of such 
liquid layers is strongly influenced by even very small 



stresses. 11 ' 12 Based on these observations, the embrit- 
tlement of Al by liquid Ga is considered to be associ- 
ated with the fast penetration of liquid gallium, result- 
ing in rapid changes in grain boundary structure/bonding 
and in many properties of the polycrystal (e.g., fracture 
strength). 12 

Thermodynamically, wetting of a grain boundary by 
a liquid metal should be expected when the spreading 
coefficients S satisfies: S — ■job — Z'JSL > 0, where 
7gs and jsl are the free energies of the grain boundary 
and solid-liquid interface, respectively. 13 However, ther- 
modynamic arguments do not explain the liquid channel 
morphology, the Ga penetration kinetics, or the atom- 
istic mechanism of Ga penetration. The anomalously 
fast, time-independent penetration rate (several /im/s at 
room temperature 4-6 ' 10-12 ) of very long nanometer-thick 
liquid films cannot be explained in terms of the classical 
Mullins grain boundary grooving 1 nor by normal grain 
boundary diffusion. 15 

Various models have been proposed to explain the ki- 
netics and atomistic mechanisms by which the liquid 
phase penetrates quickly along grain boundaries, includ- 
ing corrosive dissolution, 16 mixed diffusion-dissolution, 17 
self indentation-internal solution, 18 coherency stresses, 19 
amorphous grain boundary/liquid transformation, 20 and 
others. 21 However, even for the well-known Al-Ga LME 
couple case, many fundamental questions remain. For 
example, the degree to which Ga penetration along grain 
boundaries occurs in the complete absence of internal 
or external stresses is still unsettled. Furthermore, the 
role of stress on grain boundary penetration is not well 
understood. Further issues include whether LME is es- 
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sentially "replacement-like" (Ga atoms replace Al atoms 
at the grain boundary and the Al atoms are transported 
away) or "invasion-like" 10 (Ga atoms insert into the grain 
boundary without replacing Al atoms). The latter pro- 
cess must generate stresses, which could be relieved by 
plastic deformation as the liquid layer thickens. Clearly 
the mechanisms by which grain boundary penetration oc- 
curs remain unsettled. 

The penetration of a liquid phase along grain bound- 
aries is a complex phenomenon, involving several dif- 
ferent types of simultaneous processes; e.g., disso- 
lution/reprecipitation, liquid groove formation, grain 
boundary diffusion, and grain boundary segregation. Be- 
cause of the interplay between the underlying phenomena 
that occur in LME, it has been difficult to design experi- 
ments that can be easily interpreted to understand which 
processes control the phenomena and which are simply 
parasitic. We approach LME by performing molecular 
dynamics (MD) simulations of an Al bicrystal in contact 
with liquid Ga (with and without an applied stress) . The 
advantage of a simulation approach is that it is much 
easier to interrogate both microscopic and macroscopic 
events that occur during LME with a resolution rarely 
accessible to experiments (Angstroms and picoseconds). 
As a first step, we tune a set of embedded-atom method 
(EAM) potentials to reproduce the experimental solid- 
liquid Al-Ga binary phase diagram. We investigate how 
Ga penetrates along the grain boundaries, the degree to 
which Al dissolves into the liquid Ga, the relative dis- 
placement of the grains, the stress distribution within 
the solid, and physical parameters that are useful for in- 
terpreting the penetration mechanism (e.g., grain bound- 
ary diffusivity, segregation of Ga to the grain boundary, 
elastic constants). Based on these simulation results, we 
propose a new mechanism for LME and compare it with 
general trends gleaned from a series of experimental stud- 
ies in the literature. A brief report on this topic, focused 
upon different grain boundaries, was reported in Ref. 22. 



II. 



SIMULATION METHOD 



A. Interatomic potentials for Al and Ga 

A first step in the simulation of LME in the Al-Ga sys- 
tem is the development of a reasonable description of the 
atomic interactions. While several potentials have been 
developed for elemental Al (Ref. 23-29 and one for Ga 
(Ref. 30), we are unaware of the development of poten- 
tials that were optimized to reproduce the properties of 
the binary Al-Ga system. The Al potentials listed here 
are of the EAM form, while the Ga potential is a mod- 
ified embedded-atom method (MEAM) potential (that 
includes angular terms) 30 . While the MEAM Ga poten- 
tial reproduces many of the properties of Ga, we focus 
here on the computationally more efficient EAM class of 
potential for both Al and Ga. To this end, we develop 
Al-Ga EAM potentials of the EAM type originally de- 



veloped by Mei and Davenport 31 for Al. The parameters 
in the Mei and Davenport potentials were fitted to a se- 
ries of basic properties for which data were available, in- 
cluding the cohesive energy lattice constant, and elastic 
constants of the face-centered cubic (fee) crystal. The 
advantage of this potential is that it is analytical and 
hence easily adjustable to obtain desired thermodynamic 
properties and extendable to binary alloys. 

The total energy of the system is given by the usual 
EAM form: 
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where F Si (p) is the energy associated with embedding 
atom of type Si in a uniform electron gas of density p and 
(f) SiSj (r) is a pairwise interaction between atoms of type 
Si and Sj separated by a distance r. The pair potential 
term (f>(r) is chosen to take the form 
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where r is the first-nearest-neighbor distance in the ref- 
erence structure and <f>o, 7 and 5 are fitting parameters. 
The electron density is given by 



Pi 



(3) 



where the atomic electron density f(r) was numerically 
determined 31 so that the electron density of each atom 
in the reference structure, p re f, satisfies the following re- 
lation: 



Prof = ^2 z m f(r) = p e exp 



m— 1 



(4) 



Here, z m is the number of mth-nearest-neighbor atoms 
in the reference structure and the parameter /3 quanti- 
fies the distance over which the electron density decays 
away from an atom position. Accordingly, the embed- 
ding function takes the form 27 (including interactions up 
to third- nearest-neighbors): 



l-|ln 



F(p) = -E 
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so that the total energy of reference system func- 
tion of dilation satisfies the universal binding energy re- 
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TABLE I: Basic input data used in fitting the parameters for 
Ai[27] and Ga[30]. The numbers in parentheses are caiculated 
values. 



properties 


Al 


fcc-Ga 


E c (eV/atom) 


3.39 (3.39) 


2.897 (2.897) 


a (A) 


4.05 (4.05) 


4.247 (4.247) 


B (GPa) 


76 (76) 


52 (52) 



lation 



32 
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exp 
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-a | 1 



(6) 

with a — ^J%BVl/ Eq : where Eo, B, and fl are the co- 
hesive energy, bulk modulus, and atomic volume, respec- 
tively. The potential interactions were truncated between 
the third- and forth-nearest-neighbor shells of a static fee 
crystal. 

For Al, we used the parameter set suggested by Stur- 
geon and Laird 33 (which is a modification of those origi- 
nally obtained by Mei and Davenport 27 to reproduce the 
appropriate melting point of Al). Based on the same for- 
malism as Mei and Davenport, we fitted EAM potential 
parameters for Ga using data obtained from first prin- 
ciples with a reference fee lattice 30 (see Table I). Pa- 
rameters Eq, ro, and a were directly determined from 
the cohesive energy E c , lattice parameter a, and bulk 
modulus B, while 7 and 5 were obtained by optimization 
against the other elastic constants, vacancy formation en- 
ergy, and melting point. The potential parameters for Al 
and Ga are shown in Table II. 

Unfortunately, the new potential for Ga predicts that 
the fee crystal structure has a lower energy than the ex- 
perimentally observed All a-Ga structure. This is not 
a particularly important problem, however, since in the 
present situation, we focus only on liquid Ga. Fortu- 
nately, the liquid properties of our model Ga potential are 
in good agreement with experiment. The melting point 
of fee Ga (obtained via microcanonical ensemble molec- 
ular dynamics simulations of solid-liquid coexistence 34 ) 
is 305 K which is only about 3% higher than the exper- 
imental value for a-Ga. Figure 1 shows the structure 
factor S(k) and the pair correlation function g(r) at 959 
K obtained using this potential in molecular dynamics 
simulation and from experiment. 35 Except for the first 
peak, the calculated S(k) is in good agreement with ex- 
perimental data both in peak heights and positions. The 
pair correlation function shows a broad first peak cen- 
tered at ~ 2.8 A and a very weak second peak at ~ 5.4 
A, in good agreement with the experimental data. How- 
ever, the first peak obtained using the present potential 
is slightly broader than that from experiment. This small 
difference is consistent with that found using the MEAM 
potential. 30 The diffusivity of liquid Ga (Fig. 2) at high 
temperature obtained using the present potential is in 
better agreement with experiment 36 than that obtained 
using the MEAM potential. 
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FIG. 1: (a) Structure factor and (b) radial distribution func- 
tion for liquid Ga at 959 K for the present EAM potential 
and experiment. 35 
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FIG. 2: Arrhenius plot of the self-diffusivity of liquid Ga ob- 
tained using the present EAM potenial, the MEAM potential 
by Baskes et al., 30 and experiment. 36 



B. Potentials for the Al-Ga alloy system 

The simple analytic, EAM potential form can be easily 
extended to multi-component systems. For EAM binary 
alloys, we must fix seven functions fAi(r), fca^r), Fai(t), 
Fca(r), <pAi-Ai(r), (f>Ga-Ga(r) and (f> A i-Ga(r)- The first 
six of these are transferable from the two monatomic sys- 
tems, Al and Ga. 

We can determine the remaining function, (j>Ai-Ga(r), 
by fitting parameters O , r , 7, and S to alloy proper- 
ties. Here, we choose this function to reproduce the bi- 
nary solid-liquid alloy phase diagram of the Al-Ga sys- 
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TABLE II: Parameters for the Al-Ga EAM potential. Parameters for Al were developed by Sturgeon and Laird. Eo and 4>o 
are in units of eV, ro is in units of A, and the other parameters are dimensionless. 





Eo 


<j>0 


ro 


a 


/3 


7 


8 


Al 


3.39 


0.1318 


2.8638 


4.60 


7.10 


7.34759 


8.45 


Ga 


2.897 


0.064 


3.003 


4.42 


7.10 


7.8 


5.2 


Al-Ga 




0.075 


2.933 






5.7 


6.3 
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Experiment 
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FIG. 3: Calculated and experimental 40 phase diagram of Al- 
Ga binary alloy. 
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FIG. 4: Arrhenius plot of Al diffusivity in molten Ga for the 
present EAM alloy potential and experiment. 41 



tern. Of particular importance in the LME of Al by Ga, 
is the solubility of Ga in Al and vice versa. Our earlier 
work on general trends in binary phase diagrams as a 
function of atomic interactions 37 provides some guidance 
here. The equilibrium solid-liquid phase boundaries are 
monotonic functions of the fitting parameters <po, ro, 7, 
and 6. We initially chose 4>o to be the geometric mean of 
the monatomic pair potentials and the parameters ro, 7, 
and S to be the arithmetic mean of those of Al and Ga and 
then made small adjustments to them to obtain reason- 
able agreement with experiment. The solid-liquid alloy 
phase diagram was obtained using the Gibbs-Duhem in- 
tegration method, as described in Ref. 38 and 39. The po- 
tential parameters are given in Table II. Figure 3 shows 
the calculated and experimental 40 phase diagram for the 
Al-Ga alloy system. The agreement between the two is 
very good except near the pure Ga end of the phase di- 
agram(where it is still pretty good). 

We also measure the diffusivity of Al in liquid Ga using 
the present alloy potential. As seen in Fig. 4, the pre- 
dicted diffusivity has the same activation energy as that 
found in experiment and the absolute diffusivities agree 
to within a factor of two. This is considered excellent 
agreement for diffusivity measurements. 



C. Molecular dynamics simulations of LME 

All of the simulations described below were performed 
in a three dimensional Al bicrystal sample that is in con- 
tact with liquid Ga. Figure 5 shows a schematic illustra- 
tion of the simulation cell and initial atomic configuration 



of the grain boundary. We first constructed a periodic Al 
bicrystal by joining two unrelaxed and rotated fee crys- 
tals and a separate liquid Ga phase (the structure was 
obtained by heating to above its melting temperature). 
After relaxing these structures at the desired tempera- 
ture by MD, we brought the two systems together into a 
single simulation box with a liquid layer on top of the 
bicrystal in such a way that the grain boundaries in- 
tersect the solid-liquid interface. We imposed periodic 
boundary condition in x (normal to grain boundaries) 
and y directions, fixed one or two bottom layers of atoms 
in the solid to maintain the bicrystal structure and left 
the top surface free (i.e., there is a vacuum above the liq- 
uid). Because of the periodic boundary conditions, two 
identical grain boundaries (noted as 'GB I' and 'GB IF 
in Fig. 5) are contained in a simulation cell and the solid 
bicrystal is effectively equivalent to a bamboo-like poly- 
crystal with periodic grain boundaries separated by grain 
size dcB in x direction. 

Since the liquid Ga penetrates deep into the grain 
boundaries during the simulations (~ 70 ns) , we employ a 
relative thick bicrystal (at least 40 nm). Furthermore, in 
order to minimize interactions between the two adjacent 
liquid grooves, the simulation cell dimension in the x di- 
rection was also large (~ 66 nm) . Finally, the simulation 
cell dimension in the y direction was chosen to be only 5 
cubic unit cell lattice parameters (~ 2 nm - greater than 
the range of the EAM potentials). A typical simulation 
cell contained ~ 350, 000 atoms. 

Our simulation study focuses primarily on the 
£5 36.9°/ [010] symmetric tilt grain boundary. We chose 
this special boundary as a starting point because it is 



5 




FIG. 5: Simulation cell containing two Al grains in contact 
with liquid Ga. The enlarged picture shows the atomic struc- 
ture of the £5 36.9 °/[010] symmetric tilt boundary. 
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FIG. 6: Stress versus strain for a periodic bicrystal containing 
a £5 36.9°/[010] symmetric tilt boundary at 600 K. 



particularly well-studied and has a simple structure. 

The simulations were performed under uniform, tensile 
loading conditions. In a macroscopic sample, most of the 
load is carried by the solid, far from the surface. However, 
in atomistic simulations, the sample is necessarily small 
and there is no material far from the surface; hence, the 
effect of boundary conditions for constant stress (NPT) 
and constant strain (NVT) simulations can be quite dif- 
ferent. Since the goal is to model macroscopic loading 
conditions, we performed the simulations under constant 
strain conditions, where the strain was chosen to provide 
the desired uniaxial stress in the bulk bicrystal. To ac- 
complish this, we applied a constant strain by hxing the 
length of the simulation cell in the x and y directions. 
Note that the elastic modulus of solid bicrystals depends 
on the crystalline orientation of the two grains as well 
as the grain boundary characteristics. To find the stress 
that corresponds to the fixed length simulation cell, we 
measured the stress-strain curve in a bicrystal cell that 
is periodic in all directions and contains no liquid. This 
is shown in Fig. 6. For the £5 36.9 °/[010] symmetric tilt 
grain boundary, the stress-strain relation of the bicrystal 
cells shows elastic behavior up to strain of ~ 2%. 

The MD simulations were performed in the NVT 



(canonical) ensemble using the LAMMPS (Large- 
scale Atomic/Molecular Massively Parallel Simulator) 
code. 42 ' 43 The equations of motion were integrated us- 
ing the velocity Verlet method. The simulations were 
performed at T = 600 K and with an applied uniaxial 
stress in the to 500 MPa range. Total simulation time 
was at least 70 ns (~ 3 x 10 7 At, where At =2.5 fs rep- 
resents the time step for the integration of the equation 
of motion) . All simulations in this study were performed 
on an IBM BlueGene/L machine at Princeton University 
with up to 2048 processors. 



III. PENETRATION OF GALLIUM INTO £5 
SYMMETRIC TILT BOUNDARIES WITH AND 
WITHOUT APPLIED STRESS 

To investigate the wetting of a grain boundary in Al 
in contact with liquid Ga under the influence of external 
stresses, we performed a series of MD simulations using 
the geometry of Fig. 5 at T = 600 K. Figure 7 (a) shows 
the Ga penetration along the £5 36.9°/[010] symmetric 
tilt boundaries in the absence of an applied strain. At 
the beginning of the simulation (within less then 10 ns), 
some Al atoms near the grain boundaries dissolve into the 
liquid Ga and liquid grooves form at the intersections of 
the grain boundaries and the solid- liquid interfaces. Note 
that at T = 600 K, the solubility of Al in liquid Ga pre- 
dicted by the phase diagram is ~ 30 mole percent. The 
amount of Al that can dissolve into the liquid is limited 
by the finite quantity of liquid Ga present in the system. 
As the liquid phase approaches saturation, the shape of 
the solid-liquid interface evolves very slowly and the rate 
of growth of the liquid groove decreases. However, below 
the root of the liquid groove (i.e., where the liquid meets 
the grain boundary), Ga atoms continue to penetrate into 
the grain boundaries. Examination of the simulation cell 
at the atomic level shows that the Ga distribution ap- 
pears to be diffusion-like rather than appearing as a liq- 
uid that abruptly terminates at a crack tip. This is Ga 
penetration along the grain boundary. Overall morpho- 
logical and composition equilibrium is not yet obtained 
during the course of our 70 ns simulations. 

Figures 7(b) and (c) show the formation of liquid 
grooves and Ga penetration at constant strains of about 
0.65 % and 1.3 %, respectively. The imposed strains 
correspond to uniaxial stresses of a xx « 250 MPa and 
(7 xx ~ 500 MPa (i.e., normal to the nominal grain bound- 
ary plane), respectively, as described in section II C. Al- 
though the liquid groove shape and wetting angle are 
nearly the same in Figs. 7(a)-(c), Ga penetration into 
the grain boundaries at the root of the liquid groove was 
greatly promoted by the application of the strain. 

We also performed similar simulations for several other 
types of grain boundaries and found that the Ga pene- 
tration behavior and the effect of applied stresses can be 
very sensitive to grain boundary crystallography. For ex- 
ample, we were unable to see any grain boundary wetting 
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FIG. 7: Atomic scale images of liquid metal penetration in an Al bicrystal in contact with liquid Ga at t — 10, 30, and 50 ns 
(from top to bottom) at constant strains corresponding to applied stresses <j xx of (a) 0, (b) 250, and (c) 500 MPa. The atoms 
shown in grey represent Al atoms and those in black are Ga. 



within our 50 ns simulations in low angle grain bound- 
aries or the £5 36.9 °/[100] symmetric twist grain bound- 
ary. Since the £5 36.9°/[010] symmetric tilt boundary 
is a well-studied grain boundary in the literature and 
exhibits remarkable rate of Ga penetration in the simu- 
lations, we limit our focus to this grain boundary in the 
remainder of this paper. 

Figure 8 shows the effect of applied stress on Ga pen- 
etration rates. In this figure, we plot the product of the 
grain boundary width 5 and the Ga concentration Xq^ 
(Ga atoms per volume) versus distance along the sym- 
metric tilt boundaries. (The Ga concentration is mea- 
sured in thin slices through the sample that are perpen- 
dicular to the nominal grain boundary.) At the initial 
stage of liquid groove formation (t < 5 nm) , there is lit- 
tle effect of applied stress on the penetration profile. This 
regime is dominated by dissolution. We can arbitrarily 
define the Ga penetration depth by noting the depth at 
which the Ga concentration exceeds a fixed value at each 



time. We mark the depth at which the grain boundary 
Ga concentration SXq^ exceeds about half a monolayer 
(~ 6 atoms/nm 2 ) at each time in Fig. 8. As the grooves 
deepens, stress facilitates the rate of grain boundary pen- 
etration. For example, Ga penetrates at least 10 nm 
along the grain boundaries between t = 10 and 50 ns 
when an applied stress is present, but less than 5 nm 
when no stress is applied. 

We plot the Ga penetration depth L versus time t in 
Fig. 9. In the absence of an applied stress, the rate at 
which Ga penetrates down the grain boundary (slope in 
Fig. 9) decreases with time. However, when a stress is 
applied, the rate of Ga penetration appears to be nearly 
independent of time (the data in Fig. 9 fall on straight 
lines). [In the lower applied stress (<~ 250 MPa) case, only 
one grain boundary (GB I) shows a time-independent Ga 
penetration rate. This is likely attributable to the fact 
that for the present very small grain size, Ga penetration 
along one grain boundary relieves the stress at the other 
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Depth (nm) 

FIG. 8: Ga penetration (8Xq^) profiles along the symmetric 
tilt grain boundaries with different applied stresses at (a) t = 
10 ns, (b) t = 30 ns, and (c) t — 50 ns. The solid lines are 
cubic spline-fits and the vertical bars provide a measure of 
the Ga penetration depth, arbitrarily denned as the depth at 
which the Ga concentration is 6 atoms/nm 2 (solid horizontal 
line). 
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FIG. 9: Ga penetration depth versus time with different ap- 
plied stresses. This depth is the distance from the surface 
where the Ga concentration SXq^ = 6 atoms/nm 2 (see the 
horizontal line in Fig. 8). The two sets of data at 250 MPa 
correspond to two different grain boundaries in the same sim- 
ulation. 



grain boundary.] Clearly, stress changes the fundamental 
nature of Ga penetration down grain boundaries in Al. 
The time independence of the Ga penetration rate shows 
that the Ga is not simply undergoing random walk dif- 
fusion on the grain boundary (L oc t 1 / 2 ) nor by normal 
grain boundary grooving 15 (L oc t 1 ^ 3 or L oc t 1 / 4 for bulk 
or surface diffusion control, respectively), but is strongly 
driven [L oc t). The Ga penetration rate increases with 
the magnitude of the applied tensile stress and increasing 
grain size doB- 
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FIG. 10: Fraction of Al in the liquid versus time for simula- 
tions performed under different applied stresses. 



IV. EFFECTS OF STRESS ON GRAIN 
BOUNDARY PROPERTIES 



To investigate the origin of the stress effect, we ana- 
lyzed several, potentially relevant, physical properties as- 
sociated with the grain boundary as a function of applied 
stress. Figure 10 shows the time evolution of the mole 
fraction of Al in the liquid phase in the MD simulations 
of Ga penetration reported above. This quantity is re- 
lated to the rate of Al dissolution (primarily) at the grain 
boundary. The mole fraction of Al in the liquid increases 
with time from near zero, but appears to saturate at late 
times during the MD simulations. The equilibrium solu- 
bility of Al in liquid Ga predicted by the phase diagram 
calculation should be ~ 30 mole percent. This is very 
close to that observed by the end of the 70 ns simulations. 
However, the three curves in Fig. 10 are nearly indistin- 
guishable from each other. This demonstrates that stress 
has little effect on the rate or magnitude of Al dissolu- 
tion. This may not be surprising since the stress (normal 
traction) at the solid-liquid interface should be zero even 
in the presence of an applied stress. Furthermore, appli- 
cation of a 500 MPa uniaxial stress should only lower the 
melting point by ~ 1° and increase the solubility of Ga 
in Al by much less than 1% (based on the increase in the 
energy of the solid phase with strain and on the phase di- 
agram). Therefore, we cannot attribute the stress-effect 
to stress-driven changes in the solubility of Al in liquid 
Ga. 

Another possible explanation for the role of an applied 
stress on the rate of Ga penetration down grain bound- 
aries in Al may be associated with the effect of stress on 
Ga diffusion along the grain boundaries in Al. Unfortu- 
nately, it is difficult to measure the diffusivity of Ga along 
an Al grain boundary because the number of Ga tracer 
atoms would be too small to achieve sufficient statistical 
accuracy at low Ga concentration and, at high concen- 
trations, the grain boundary pulls apart in the presence 
of an applied stress (at 600 K). If the impurity diffusivity 
is correlated with the self-diffusivity via a vacancy mech- 
anism, 44 the Al self-diffusivity along the grain bound- 
ary is relevant. Figure 11(a) shows the in-plane mean- 
squared displacement (Ar(i) 2 )(= (Ax(t) 2 ) + (Ay(t) 2 )) 
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FIG. 12: (a) An atomic scale image Irom a semi-grand canoni- 
cal Monte-Carlo simulation of a £5 36.9 °/[010] symmetric tilt 
boundary in the absence of an applied stress at 600 K. The 
white and black circles represent Al and Ga atoms, respec- 
tively, (b) The Gibbsian grain boundary excess Tea of Ga 
versus the applied stress. 



FIG. 11: Grain boundary thickness - self-diffusivity product 
for the E5 36.9°/[010] symmetric tilt boundary at 600 K. 
(a) The mean-square displacement of the Al atoms in the 
grain boundary at different uniaxial stresses as a function of 
time, (b) The logarithm of the grain boundary thickness - 
self-diffusion coefficient as a function of the applied stress. 



for Al atoms in the X5 symmetric tilt boundary as a func- 
tion of time at different uniaxial stresses, u X x, in an ef- 
fectively infinite Al bicrystal (i.e., no Ga/no liquid). The 
mean squared displacement is a linear function of time. 
Figure 11(b) shows the the boundary width (SDqb) - 
self-diffusion coefficient product obtained from the slopes 
of the curves in Fig. 11(a). These results suggests that 
the grain boundary self-diffusivity is very nearly indepen- 
dent of applied stress (to within the accuracy of the sim- 
ulations). Therefore, the effect of stress on grain bound- 
ary diffusion is an unlikely source for the difference in Ga 
penetration rates observed with different applied stresses. 

Next, we investigated the effect of stress on the seg- 
regation of Ga to the grain boundary. To this end, we 
perform a series of Monte Carlo simulations of a bicrys- 
tal (no liquid) in a semi-grand canonical ensemble where 
we fix the chemical potentials of Al and Ga such that 
the bulk concentration of Ga is ~ 10% at T = 600 
K. To simulate segregation in the presence of a uniax- 
ial stress, we adopted the isotension-isothermal ensem- 
ble. 45 Figure 12(a) shows the atomic structure of an Al- 
Ga bicrystal with symmetric tilt grain boundaries in the 
absence of an applied stress. Clearly, Ga segregation to 
the grain boundary does occur and leads to a relatively 
disordered grain boundary structure at T = 600 K. We 
measure the Gibbsian grain boundary excess Tc a of Ga 
at the grain boundary from the concentration profile for 



several applied stress, as shown in Fig. 12(b). The grain 
boundary excess is the quantity of solute present per unit 
area of flat interface in excess of the quantity of solute 
that would be present if there were no grain boundary 
(i.e., the bulk concentration). The grain boundary may 
be defined as T Ga = (N tota i - C bu i k V to tai)/A GB: where 
Ntotai is the total number of solute atoms in the cell vol- 
ume Vtotai, Cbuik is the bulk concentration of solute, and 
Aqb is the grain boundary area in the volume. As shown 
in Fig. 12(b), applied stresses have little impact on the 
degree to which Ga segregates to grain boundaries in Al. 



V. DISTRIBUTION OF STRESSES 

Since stress has little effect on many of the basic physi- 
cal parameters that could, in principle, be responsible for 
the stress effect in LME, we examine the distribution of 
the stress within our Al-Ga bicrystal system. Figure 13 
shows the time evolution of the stress distribution for 
a xx within the system at constant strains of 0, 0.65% 
(a xx « 250 MPa), and 1.3% {a xx « 500 MPa). Fig- 
ure 13(a) shows that in the absence of an applied strain, 
the stresses in the system are small and random. How- 
ever, when a strain is applied, we note the formation of 
one [Fig. 13(b)] or more [Fig. 13(c)] stress concentration 
patterns at the grain boundary. These patterns consist of 
a dark (large compressive) region above a light (large ten- 
sile) region. This stress patterns are suggestive of stress 
patterns expected for edge dislocations with a Burgers 
vector perpendicular to the boundary plane. 

We can use linear elastic theory to predict the stress 
field associated with such an edge dislocation in a linear 
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(a) a xx = MPa (b) a xx = 250 MPa (c) a xx = 500 MPa (d) Dislocation model 

FIG. 13: Stress distributions (a xx ) at t = 10, 30, and 50 ns (from top to bottom) for simulations performed at T = 600 K at 
constant strains corresponding to applied stresses a xx of (a) 0, (b) 250, and (c) 500 MPa. (d) The stress distribution predicted 
from the dislocation model (see the text) is shown for comparison with (c) . 



elastic half-space using the approach of Head 46 . This 
stress field a xx corresponding to a set of edge dislocations 
with Burgers vector b parallel to the free surface z=0, 



located at points in the xz-plane (see Fig. 5) is 

given by 
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where E and v are the Young's modulus and Poisson's ra- 
tio, respectively. Figure 13(d) shows the predicted stress 
field corresponding to Fig. 13(c), where we have assumed 
that the positions of the dislocations in Fig. 13(d) are 
the same as those observed in Fig. 13(c). In making this 
comparison, we have determined the magnitude of the 
Burgers vector in Fig. 13(d) by adjusting it to obtain best 
correspondence with the stress distributions in Fig. 13(c). 
We find that the best fit is obtained for |6| ~ 0.5 A, which 
is much smaller than the Burgers vector of a lattice dis- 
location, ~ 2.86 A. 

No dislocations were formed in the case of zero applied 



stress [Fig. 13(a)] or at one of the two grain boundaries 
in the a xx = 250 MPa [Fig. 13(b)]. It is interesting to 
note that it was only in exactly these two cases that the 
Ga penetration depth versus time plot (Fig. 9) were sub- 
linear (i.e., the penetration rate decreases with time). 

Examination of Fig. 13(b) shows that a dislocation, 
once formed, "climbs" down along the grain boundary 
at a nearly constant rate. Examination of Fig. 13(c) 
(larger strain) shows that the dislocations "climb" down 
the grain boundaries at the same constant rate as the sin- 
gle dislocation in Fig. 13(b) (low strain). However, in this 
case, once the first dislocation has moved some distance 
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from its point of origin (below the liquid groove root), 
second and third dislocations are nucleated one after an- 
other (boundary on the left in Fig. 13(c)), and "climb" 
down the grain boundaries too, leading to equally spaced 
dislocations that all move at the same rate. Therefore, 
we can conclude that these special grain boundary dis- 
locations only form above a critical applied strain/stress 
and "climb" down the grain boundaries at a constant 
rate that is independent of the magnitude of this strain. 
Increasing applied strain simply results in the formation 
of more dislocations with shorter incubation time. 

In order to understand the relaxation of the stresses, 
we also analyzed the atomic-level displacement fields. 
Figure 14 shows the Ga concentration profiles at t = 50 
ns (contour plots) and the displacements in the solid that 
occurred between t — 10 and 50 ns (i.e., the displacement 
vectors are measured as the atom positions at 10 ns to the 
positions of the same atoms at 50 ns) at applied strains of 
0, 0.65% (a xx w 250 MPa), and 1.3% {a xx w 500 MPa). 

In the absence of an applied stress, only small displace- 
ments were observed and these displacements were dis- 
tributed randomly over the entire solid (consistent with 
the stress distributions shown in Fig. 13(a)). When a 
stress is applied, the atomic displacements are consid- 
erably larger than in the absence of an applied stress 
and a distinctive displacement vector pattern emerges 
[see Figs. 14(b) and (c)]. The atomic displacements near 
the grain boundary groove tip are directed away from 
the grain boundary. This is likely associated with Ga 
atoms being inserted into the grain boundary to reduce 
the tensile stresses there. The largest atomic displace- 
ments observed within the grains adjacent to where Ga 
has thoroughly penetrated the nearby grain boundary. 

Additional simulations, not shown here, demonstrate 
that doubling the grain size doubles the magnitude of 
the displacements. This can be traced to the fact that 
the present simulations were performed under fixed-grip 
loading (to a particular bulk stress value). This implies 
that doubling the grain size doubles the strain energy 
stored in the sample. To relieve the same stress in a 
system with twice the grain size requires twice the grain 
boundary opening displacement. This is consistent with 
the observation that increasing the grain size leads to 
a dramatic increase in the quantity of Ga at the grain 
boundary. 47 

We also analyzed the relative displacement between a 
pair of points in adjacent grains as an indicator of grain 
boundary opening; the relative displacement A5 a b(t)(— 
S a b(t) — S a b(0)) is the change of the distance 5 a b{t) be- 
tween the two points indicated in Fig. 5 at time t with 
respect to the distance at zero time. This displacement is 
a measure of the grain boundary opening distance. This 
type of analysis was used earlier in experimental observa- 
tions of LME (e.g., see the Al-Ga experiments in Ref. 10). 
Figure 15 shows the relative displacement AS a b(t) of a 
pair of points from adjacent grains and the effective Ga 
layer thickness WQa as a function of time. (The effective 
thickness of the Ga layer quantifies the quantity of Ga at 



the grain boundary and is defined as WQa — NQ a Q/AcB, 
where Nc a is the number of Ga atoms in a slab perpen- 
dicular to the grain boundary times the atomic volume 
of Ga and divided by the cross-sectional area of the 
slab Aqb-) AS a b and wc a are shown as a function of 
time in Fig. 15, as measured at z = 15 nm below the 
initial surface. Clearly, the relative displacement S a b is 
initially zero then, at a finite time, gradually increases 
as Ga penetrates into the grain boundary. The effective 
Ga layer thickness wc a is also initially zero, but then 
increases abruptly at a time later than where the rela- 
tive displacement begins to grow. The abrupt increase 
in the effective Ga layer thickness occurs when a dis- 
location passes the measurement depth z = 15 nm, as 
observed in Fig. 13(c). Subsequent abrupt rises in wc a 
correspond to the passage of additional dislocations. In- 
terestingly, the effective Ga layer thickness is comparable 
to the grain boundary opening distance. A similar evolu- 
tion of A5 a b(t) and WGa{t) occurs at points farther from 
the surface (not shown) at later times (set by the disloca- 
tion climb velocity). Although the time and length scales 
are different, the shape of these curves from the simula- 
tions in Fig. 15 are very similar to those measured ex- 
perimentally by Ref. 10. The following picture emerges. 
Ga diffuses down the grain boundary leading to some 
grain boundary opening. At a critical opening, a dis- 
location forms and starts to climb downward, trailing a 
Ga concentration tail behind it. This Ga concentration 
at the boundary and the dislocation itself lead to some 
grain boundary opening ahead of the dislocation (and Ga 
concentration tail). The grain boundary opening, dislo- 
cation and Ga tail move down the grain boundary at the 
same rate. In the next section, we propose a mechanism 
by which all of this (i.e., LME) occurs. 



VI. DISCUSSION 

We performed a series of MD simulations of an Al 
bicrystal in contact with liquid Ga and examined the pen- 
etration of Ga along the Al grain boundaries in the pres- 
ence and absence of an applied stress. Even though the 
present simulations were relatively large, both in terms 
of the number of atoms modelled and simulation time, 
our atomistic simulation approach to LME, nonetheless, 
has some limitations. For example, the simulations were 
only able to describe the first ~ 70 ns of the LME pro- 
cess and the extent of the liquid metal penetration was 
limited to <~ 40 nm. This should be compared with many 
experiments in which thin (nanometer-thick) liquid films 
penetrate up to several hundred micrometers over time 
scales ranging from seconds to hours. Because of these 
limitations, the stresses applied in the present simulation 
had to be much larger than those applied in experiment 
in order to observe LME. The relatively small size of our 
simulation cells also precluded observations of significant 
plastic deformation (dislocation nucleation in small vol- 
umes is difficult and there were no pre-existing disloca- 
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FIG. 14: Ga concentration profiles (contour plots of mole fraction Xa a ) and displacement fields, measured between t= 10 
and 50ns, for the bicrystal with an applied stress of (a) 0, (b) 250, and (c) 500 MPa. The displacement vectors correspond 
to averages over several hundred atoms and are shown magnified by a factor of ten for better resolution. In order to limit 
consideration to displacements associated with elastic deformation of the solid, we excluded atoms for which the displacements 
exceeded 5 A (i.e., primarily diffusive hops) and, as a result, no arrows are plotted in the grain boundary region. 
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FIG. 15: Relative displacement between two points in ad- 
jacent grains versus time with an applied stress of 500 MPa. 
The distance between a pair of points was measured at z = 15 
nm below the initial surface (see Fig. 5). The effective Ga- 
rich layer thickness WGa versus time measured at the same 
position is also shown. 



tions), even though plastic deformation would occur in 
macroscopic samples at this stress level. However, the 
fact that significant plastic deformation does not occur 
during our simulations was fortuitous since experiments 
show that LME occurs in Al-Ga at stresses which are too 
small to cause macroscopic deformation of the Al. 

The molecular dynamics simulations can complement 
continuum models and experiments. First, simulations 
can provide an atomic scale view of the dynamics of the 
system, grain boundary structure, and stresses and com- 
position near the advancing liquid film layer front. Such 
fundamental, high resolution data are rarely accessible in 
experiments: for example, although the synchrotron ra- 
diation X-ray microradiography 10-12, 48-51 is very useful 
for in-situ observation of sub-micrometer Ga wetting lay- 
ers, the very early stages of the wetting process involves 
a Ga layers that are only a few monolayers thick and are 
not visible with this technique (typically, a layer thick- 
ness of at least ~ 5 nm is required for detection). The 



molecular dynamics simulations provide a unique tool for 
understanding which materials properties are important 
in LME and for identifying which processes are rate- 
controlling. This is accomplished by careful design of 
idealized set of rigidly maintained conditions to empha- 
size certain effects and exclude competing factors which 
may influence the LME. While no interatomic potential 
is perfect, the EAM potentials for the Al-Ga alloys were 
shown to successfully reproduce many of the important 
solid-liquid properties in this alloy system. To comple- 
ment the known behavior of these potentials, we also 
performed an analysis of several other properties of this 
system that may affect LME, including grain boundary 
diffusivity and the tendency for Ga grain boundary segre- 
gation. This gave us the ability to exclude certain models 
for LME from further consideration. 

Our simulation results may be directly compared with 
many experimental observations in the literature. This 
is, in part, associated with the fact that Al-Ga is one 
of the most widely experimentally studied LME systems 
(both polycrystals and bicrystals). Among these, the re- 
cent series of TEM studies by Hugo and Hoagland, 4-6 and 
synchrotron radiation micro-radiography studies by Lud- 
wig and Pereiro-Lopez et aZ. 10-12,48-51 are of particular 
interest because they are both quantitative and provide 
microscopic observations. Although these experimental 
techniques are very different, both studies report a con- 
sistent set of features. 

Although the time and length scales of the simulations 
and experiments differ, we demonstrated that our sim- 
ulations were able to capture many of the experimen- 
tally observed trends in liquid Ga penetration of grain 
boundaries in Al. Experiments showed that the applica- 
tion of tensile stresses (as small as a few MPa) to an Al 
polycrystal sample, drastically increased the liquid metal 
penetration rate. 52 Recent synchrotron radiation X-ray 
micro-radiographic experiments 51 showed that the pene- 
tration behavior of liquid Ga along two different types 
of symmetrical tilt bicrystals of Al was greatly facili- 
tated by the application of a tensile stress of only ~ 5 
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MPa. Our simulations also showed that an applied stress 
significantly increases the rate of liquid metal penetra- 
tion, although the applied stress level used in the simula- 
tions was approximately 100 times larger than that used 
in the experimental studies (<~ 500 versus <~ 5 MPa). 
As in the macroscopic experiments, 52,53 the simulations 
demonstrated that the penetration rate and thickness of 
the Ga-rich layer increases with increasing grain size 47 
(under fixed grip conditions). Since grain sizes in most 
real materials are much larger than those used in the sim- 
ulations (~ 1 mm versus <~ 33 nm, a factor of <~ 3 x 10 4 ), 
the difference in the magnitude of the applied stress nec- 
essary to facilitate the penetration rate in the simulations 
is understandable. 

It is clear that, even in the absence of an applied stress, 
many experimental samples may be subject to resid- 
ual stresses resulting from sample fabrication, process- 
ing, polishing, or gripping. Ludwig and Pereiro-Lopez 
et al. ' demonstrated that residual stress, introduced 
by sample preparation, can affect liquid metal penetra- 
tion rates in both Al bicrystals and polycrystals. In their 
experiment, susceptibility to grain boundary penetration 
was strongly influenced by the method in which the sam- 
ple was gripped. Two different gripping procedures were 
used in their experiments: "gently screwing" together 
between metal plates and "gluing" the Al sample to Cu 
supports using silver paint. The screwed samples exhib- 
ited larger liquid film propagation and thickening rates, 
as well as a larger probability that specific grain bound- 
aries (in Al bicrystals) could be penetrated at all. Al- 
though the magnitude of these effects was not accurately 
quantified, the sensitivity of grain boundary penetration 
to gripping methods suggests that the detailed loading 
state (e.g., fixed stress or fixed strain) may play an im- 
portant role in liquid film penetration. This is consistent 
with the discussion of the simulation results presented 
above. 

In many experiments, Ga penetrates into Al bicrys- 
tals as a thin layer that lengthened at nearly a fixed 
rate, 5 ' 7 ' 10,11 ' 50 whereas in experiments in polycrystalline 
Al, the rate of propagation was more jerky and irregu- 
lar. 12 Commonly, the penetration rate is higher in poly- 
crystals than in bicrystals presumably because of higher 
residual stresses and more degrees of freedom for the rel- 
ative motion of individual grains 49 The penetration rate 
was also observed to be very sensitive to grain bound- 
ary crystallography. In-situ synchrotron observations of 
high energy 150°/[110] tilt boundaries by Ludwig et al. 11 
showed the Ga penetration rate was nearly constant with 
a velocity v s=y 1-10 fim/s. In-situ TEM observations of 
several Al bicrystal foils by Hugo and Hoagland 5 also 
showed a nearly constant Ga penetration rate, v ~ 0.5- 
7 /im/s, depending on the grain boundary crystallogra- 
phy. These two experiments were performed under nom- 
inally similar conditions. In our simulations, we found 
that the Ga penetration rate was also nearly constant, 
albeit with a velocity of 0.1-0.2 m/s which is significantly 
larger than experimentally observed in Al bicrystals. The 



discrepancy is likely attributable to the fact that our 
simulation was performed at relatively high temperature 
(T = 600 K) and with relatively large stresses (to over- 
come the time scale limitations of MD simulations) . The 
in-situ experiments were performed at room temperature 
with a stress of several MPa, while the simulations were 
performed at 600 K with a stress of at least 250 MPa. 
Another possible reason for the faster penetration rates 
in the simulations as compared with the experiments is 
that the minimum observable Ga layer thickness in the 
simulations was ^0.2 nm versus at least ~ 5 nm in ex- 
periments. 

In the present simulations, we observed that the sep- 
aration between the pair of grains that meet at a grain 
boundary (measured using markers that were some dis- 
tance from the grain boundary) increased as the Ga-rich 
layer propagated down the grain boundary (Fig. 15). A 
similar analysis was also performed in a series of exper- 
iments by Ludwig et a/., 10,11 ' where the distance be- 
tween two grains in Al bicrystals was measured by image 
correlation techniques. They observed an increase of sep- 
aration between the points by tens of nanometers during 
penetration even without an applied stress (but possi- 
bly subject to residual stresses). Although the time and 
length scales are somewhat different, these observations 
agree well with our simulation results. The discrepancy 
between the magnitudes of the relative displacements (by 
two orders of magnitude) can be explained by the same 
stress and grain size effects described above. Both the 
simulations (Fig. 15) and experiments 10 show that the 
measured Ga layer thickness uiQ a is nearly zero until a 
time well beyond the time at which significant grain dis- 
placement was observed to begin. This suggests that the 
grain boundary opening below the liquid film tip is the 
results of the bond stretching associated with the stress 
field of the advancing Ga penetration front. 

Hugo and Hoagland 5 observed generation of a mov- 
ing strain field at the Ga penetration front using TEM. 
They proposed that the penetration front acts as a line 
defect with a singular strain field. Our simulations also 
revealed that the moving Ga penetration front was ac- 
companied by a dislocation with its own unique singular 
stress pattern (Fig. 13). This dislocation seems to pre- 
cede the advancing and thickening Ga penetration layer. 
It is interesting to note that in the absence of an applied 
strain, no dislocation forms [Fig. 13(a)] and the Ga pene- 
tration rate decreases with time (Fig. 9). However, when 
a strain is applied, dislocations form and climb at a fixed 
rate [Fig. 13(b) and (c)] and the Ga penetration rate is 
time independent (Fig. 9). This suggests that the con- 
stant Ga penetration rate observed in the strained solid is 
associated with the fixed rate of "climb" of dislocations. 
Here, the applied strain plays essential roles in Ga pene- 
tration: to aid the nucleation of dislocations at the grain 
boundary and to keep the grain boundary open to allow 
fast Ga transport enough to move with the dislocations. 

Our simulations demonstrate that application of a 
stress significantly promotes liquid metal penetration, re- 
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suiting in a change from a diffusive to fixed rate pen- 
etration mode. The change of penetration kinetics is 
attributed to the nucleation and constant-rate climbing 
of grain boundary dislocations. Why do the disloca- 
tions move down the boundary? The dislocation sets 
up its own stress field; in the present geometry, it is 
compressive above the dislocation line and tensile below. 
The chemical potential along the grain boundary is pro- 
portional to the grain-boundary traction 54 or a xx and, 
hence, the chemical potential along the grain boundary 
changes abruptly at the dislocation. Ga atoms in the 
grain boundary respond by jumping quickly from above 
the dislocation line to below it. This, in turn, moves the 
dislocation down, yet preserves the stress discontinuity 
(i.e., the dislocation is intact). This explains why the 
dislocation climbs down at a fixed rate. How fast does 
the dislocation climb? This can be determined by solving 
the coupled elasticity/diffusion problem. Similar prob- 
lems were addressed by Chuang 54-56 and Gao et al. 57 ~ 59 
in the context of diffusive crack or wedge growth along a 
grain boundary subjected to an applied stress. 56,58 ' 59 The 
steady-state dislocation climb velocity V in this model 
can be approximated as 56 
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where VL is the atomic volume of the species with grain 
boundary diffusivity D g i,, kT is the thermal energy, b is 
the Burgers vector, E is the Young's modulus, and v is 
Poisson's ratio. l c is a characteristic length associated 
with the jump in stress across the dislocation and should 
be of order of the dislocation core size (i.e., a few A) 
and can be found by solving the singular coupled elastic- 
ity/diffusion problem. 56 ' 58,59 Using values for Ga in Al in 
Eq. (8) yields V « 0.1 m/s, which is consistent with the 
dislocation climb velocity in the present simulations. 

The following model for the embrittlement of Al by Ga 
emerges. First, Ga diffuses down the grain boundary in 
Al below the liquid groove root and, if the quantity of the 
inserted Ga is sufficiently large, a dislocation is nucleated 
at the grain boundary with the aid of the applied stress. 
The dislocation establishes its own stress field; in the 
present geometry; it is compressive above the dislocation 
and tensile below it. The first dislocation climbs down 
by stress-enhanced Ga hoping across the dislocation core, 
leaving a tail of Ga behind). This Ga hopping leads to 
a constant dislocation climb rate that is independent of 
the remote applied stress. Once the dislocation moves 
far enough from the groove root, another dislocation is 
nucleated. It too climbs down the grain boundary at 
the same rate, resulting in a uniform spacing of climbing 
dislocations. Each dislocation further relaxes the applied 
stress until it reaches a level too small for further disloca- 
tion nucleation. At this point, however, the stress is not 
fully relieved. This residual stress is further relaxed by 
Ga diffusion down the boundary. Since Ga weakens the 



Al bonds at the grain boundary, 60-63 this Ga penetration 
leads to boundary decohesion, grain boundary opening, 
and crack formation/propagation (i.e., LME cracking). 
This crack can be filled with liquid Ga. The propagation 
of this "crack" leads to further Ga layer thickening. This 
"crack," Ga layer at the grain boundary and dislocations 
move down the grain boundary in unison. The Ga pene- 
tration rate mirrors the dislocation climb rate and hence 
is time independent. 



VII. CONCLUSION 

Although LME exhibits a diverse set of fracture char- 
acteristics, depending on the solid-liquid metal couple, 
the penetration of nanometer-thick liquid metal films 
along the grain boundary is one of the hallmarks of the 
process that has been observed in the classical LME sys- 
tems, such as Al-Ga, Cu-Bi and Ni-Bi. We have employed 
EAM potentials optimized for Al-Ga binary alloys in the 
performance of a series of MD simulations of an Al bicrys- 
tal in contact with liquid Ga in the presence and absence 
of an applied stress. Our simulations demonstrated how 
Ga penetrates along the £5 36.9°/[010] symmetric tilt 
boundary during the early stages of LME and how an 
applied stress enhances the Ga penetration. The simula- 
tions capture many of the experimentally observed trends 
in Ga penetration of grain boundaries in Al. The key 
atomistic mechanism at the tip of the advancing Ga pen- 
etration layer was identified through analysis of displace- 
ment fields and the stress distribution within the Al-Ga 
bicrystal system. The interplay of stress and penetrating 
Ga atoms leads to the nucleation of a train of dislocations 
on the grain boundary below the liquid groove root which 
climbs down the grain boundary at a nearly constant 
rate. The dislocation climb mechanism and the Ga pen- 
etration are coupled. While the dislocations do relax part 
of the applied stress, the residual stresses keep the grain 
boundary open, thereby allowing more, fast Ga trans- 
port to the penetration front (i.e., Ga layer thickening 
process). We believe that the coupled Ga transport and 
"dislocation climb" is the key to the anomalously fast, 
time-independent penetration of Ga along grain bound- 
aries in Al. The simulations explain most of the main 
features of a series of experimental studies of LME in 
Al-Ga over the past decade. 
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